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ABSTRACT 

Calcium-aluminum-rich inclusions (CAIs) are the oldest solar system solids known in primitive meteorites (chondrites). They predate 
the other components by 1-2 Myr, and likely condensed within a short time interval, close to the Sun in the gaseous protoplanetary 
disk. Their preservation must counterbalance both the sunward drift caused by gas drag and the general inward motion of the entraining 
gas. We propose that an efficient outward transport of CAIs can be achieved by advection as a result of the viscous expansion of the 
disk, provided it is initially less than 10 AU in size, which we argue is plausible from both observational and theoretical points of 
view. Gas drag would stop this outward motion within lO' yr. However, by that time, a magnetically dead zone would have developed 
as gravitational instabilities fade away, which would trap CAIs for a significant fraction of the disk lifetime because of the reduced 
advection velocities. The dead zone would also prevent outward diffusion of subsequently condensed CAIs, contributing to their 
observed narrow age range. This preservation mechanism is independent of the outward transport scenario (before the dead zone 
formation) and a natural consequence of considering the source of turbulence in accretion disks. 
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1. Introduction 

Chondrites are the most primitive meteorites and provide 
an invaluable record of early solar system processes. They con- 
sist of millimeter-sized solids set in a fine-grained matrix, the 
most abundant of which are chondrules, silicate spherules pro- 
duced by rapid melting and cooling of precursor material (e.g. 
ConnoUy & Desch 2004). Another noteworthy, albeit rarer 



component of chondrites are refractory inclusions,m particular 
submillimeter- to centimeter-sized calcium-aluminum-rich in- 
clusions (CAIs). Thermodynamic calculations suggest that they 
condensed at temperatures of ~ 1400- 1800 K from a gas of solar 
composition, although many have since melted and/or partially 
evaporated. With an age of ~4568 Myr (e.g. Bouvier et al. 2007| l, 
CAIs are the oldest known solar system solids. Intriguingly, pri- 
mary formation of CAIs — at least those of C V chondrites — 
seems restricted to an interval of a few 10"* yr (Bizzarro et al 
|2004| l, while chondrule formation apparently began 1-2 Myr 
later and lasted ~2 Myr ( Villeneuve et al.|[2009 1. Despite some 



variation in CAI mineralogy among chondrite groups, similar- 
ities indicate that all CAIs formed in a common environment, 
likely close to the Sun, and were then transported to chondrite- 
forming regions (e.g. MacPherson| 2005) l. Outward transport is 
also needed for two CAIs identified in dust returned from comet 
Wild 2 (e.g. Zolensky etal.|2006 1. 



Transport and/or production mechanisms of CAIs must have 
been efficient because their mass fraction in some chondrites — 
a few percent, e.g. |Hezel et al.| (|2008) — is comparable with 
what in situ condensation from a solar gas would have achieved. 
Moreover, these mechanisms must offset the sunward drift in- 
duced by gas drag on solids ( | Weidenschilling| 1 977| l and the gen- 



eral inward motion of accreting gas, which partially entrains the 
solids. These combined phenomena typically entail a drift into 
the Sun within a few 10^ yr for millimeter-sized bodies. 

Proposed CAI transport scenarii either involve jets such as 
the X-Wind model ( Shu et al.|2001| ), or operate solely inside the 
disk, e.g. with turbulent difiiision (e.g. Bockelee-Morvan et al. 
2002t |Cuzzi et al[||2003| |Ciesla|[2Ul^ , gravitational instabili- 
ties (GI; Boss|2004| , or meridional circulation — an hypothetical 
outward flow around the midplane (e.g. |Takeuchi & Lin|[2002[ 
Ciesla|2009 1. As originally formulated in the literature, jet mod- 
els and meridional circulation should operate over most of the 
disk's lifetime and thus would not reproduce per se the narrow 
age range of CV CAIs. Conversely, turbulent diffusion and GI 
would be efficient only in the early phases. As Cuzzi et al. ( 2003| l 
emphasize, disk evolution definitely matters. 

In this letter, we sketch a scenario (depicted in Fig. 1) of 
CAI transport and preservation in the disk: We adopt in Sect. 2 a 
small and massive initial disk (stage 'A' in Fig. 1). As described 
in Sect. 3, viscous heating enables CAI condensation sufficiently 
far from the Sun to have CAIs entrained outward by disk expan- 
sion (stage 'B'), before they are stopped by gas drag (stage 'C'). 
As GI fade away, a dead zone forms, where the magnetorota- 
tional instability (MRI) does not operate, in which CAIs are ef- 
fectively trapped for a few Myr until chondrite accretion (stage 
'D') because of reduced advection velocity. This is the subject 
of Sect. 4. The model and its limitations are discussed in Sect. 5. 



2. Disk model 

Observations of class objects have not yet unambiguously 
identified disks, suggesting that many disks may have begun by 
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B. VISCOUS EXPANSION 





C. OUTWARD TRANSPORT STOPS 



D. CAIs TRAPPED IN DEAD ZONE 



Fig. 1. Cartoon of the proposed scenario: CAIs form in an ini- 
tially compact and massive disk (A), which undergoes viscous 
expansion, entraining CAIs outward (B), building up a leaky 
reservoir that is eventually released at large heliocentric dis- 
tances (C), while a dead zone emerges and traps CAIs until chon- 
drite accretion (D). 



being very compact (IMaury et al. |2010 Hennebelle & Ciardi 
2009| and references therein). The lower range of plausible val 



ues of the precollapse cloud's rotation velocity and size actually 
allow a centrifugal radius — the maximum stellocentric distance 
where infalling material lands on the equatorial plane — less than 
-10 AU ( [Pickett et al.p997||Zhu et al |2010|). Hydrodynamical 
simulations of protostellar collapse by Tscharnuter et al.j([2009| 
result within ~ lO'* yr in the buildup of a compact disk with a 
hot inner region extending to 4 AU. Moreover, magnetic braking 
during the collapse has been shown to be able to remove angular 
momentum efficiently enough to allow small disk sizes, depend- 
ing on the relative inclination between angular momentum and 
the magnetic field threading the cloud (Hennebelle & Fromang 
[2008i [Hennebelle & Ciardi|2009] l. 

We therefore consider an early epoch where the disk size 
is a few AU, and adopt, for greater definition, the surface den- 



sity profile of the self-similar model taken from Lynden-Bell & 
|Pringle| ( fT974l l, similarly to |Cuzzi et alT] ( |2003] l: 



Y(t,R) 



Mpjt) 

iTTRRoit)^ 



(1) 



with Mo(t) and Roit) the disk mass and characteristic radius, 
respectively, at time t. This model assumes that the viscos- 
ity V - aCj/Q (with c, the isothermal sound speed, Q. the 
Keplerian angular velocity and a a dimensionless parameter) is 
time-independent and proportional to the heliocentric distance 
R, and therefore may be written as 



V = 5 X 10'" mVs /vRau, 



(2) 



with fy a dimensionless factor, and 7?au the heliocentric dis- 
tances expressed in AU. In the disk outer regions dominated 
by the irradiation temperature (here oc R^^l^, with nominal 
rin-(l AU) = 280 K), our nominal /v = 1 corresponds to a con- 
stant a= 10"^. Rq increases linearly with time and satisfies 



t = 



RpjtY 
MRoit))' 



(3) 



the origin of times being a mathematical singularity where all 
mass resides at the center. Angular momentum conservation im- 

1/2 1/2 

plies that M^T?^' is constant. We take a nominal M^Rj^ - 
O.SMq.AU'^^, areasonable value because, in evolved phases, the 
disk, whose radius is > 10^ AU, may be expected to have a mass 



larger than that of the 
0.01-0.02 Mo ([Hayashi 



'minimum mass solar nebula" (MMSN), 
Moreover, one-zone models of 
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embedded accretion disks by |Kratter etal., (2008) yield a disk- 
to-star mass ratio of a few tenths during the infall phase. 

It is important to note, based on Eq. (|3]l, that within ~1 Myr 
the disk, though initially small, will attain 10^-10^ AU as ob- 
served for class II objects. 

3. Outward transport via viscous expansion 

In the self-similar model, the vertically averaged radial gas 
velocity is ( Lynden-Bell & Pringlej 1974) 1 
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(4) 



where the second equality uses the surface density profile given 
by Eq. (fill. Beyond the turnover radius /^switch - ^d/2, is 
outward (positive) . 

In the inner disk, viscous dissipation generates a contribu- 
tion to the local temperature that dominates that of irradiation 
by the Sun. Using |Chambers| ( [200 9 )'s formulae, but sticking to 
the prescribed viscosity profile given by Eq. (|2]i, the former is 



lO^K /, 
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with K the (constant) opacity. Initially, as shown in Fig. 2, the 
"CAI factory" ( Cuzzi et al.|2003 i extends beyond ^switch, allow- 
ing CAIs to be advected outward from where they condensed 
(stage 'B' in Fig. 1) until /^switch exits the CAI factory on a 
timescale of 10"* - 10^ yr (see Fig. 2). This alone would account 
for the narrow age range of CV CAIs. 

The CAIs cannot be entrained arbitrarily far, however. The 
gas drag in the Epstein regime indeed induces a drift with respect 
to the gas (e.g. [Weidenschilling| 1 977\ : 



^dras ~ 



7T dP PsO 

8dRp^/ 



(6) 



with P, p the gas pressure and density, and ps, a the CAI inter- 
nal density and radius, respectively. Vdiag is inward and roughly 
inversely proportional to p. Thus, for a given CAI size, at suf- 
ficiently large heliocentric distances, Vdrag dominates Vg and the 
net CAI velocities are directed inward. Therefore, CAI velocities 
in the early phase are positive over a finite range of heliocen- 
tric distances only, bounded by two "stagnation" points (i.e. he- 
liocentric distances where the CAI velocity is zero). After some 
expansion, these stagnation points merge and the velocities are 
directed inward throughout the disk from that time onward. If, in 
the self-similar model, the position of the outer stagnation point 
is taken to be stationary in time, this yields the maximum trans- 
port range of solids R^^^ and the corresponding time fmajQ 
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For fiducial disk parameters, solids can thus be transported out- 
ward for tens of AUs within ~ 10^ yr. Because the CAIs must 



' In this calculation we approximate the temperature with its irradia- 
tion contribution, as is appropriate for large heliocentric distances. 
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Lagrangian trajectories in a self-similar disk 
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Fig. 2. Space-time diagram of the self-similar model. Dashed 
lines mark the trajectory of CAIs with p^a = 1 kg/m^ and the 
black continuous line indicates the locus of /^switch- The MRI- 
inactive areas are shaded in blue (dead zone), with the 1500 K 
isotherm (a typical CAI condensation temperature) drawn as a 
continuous blue line. Regions prone to GI in the absence of vis- 
cous dissipation are shaded in red. Once the (Gl-free) dead zone 
reaches a significant fraction of the disk mass, the disk model is 
no longer self-consistent, and disk evolution is expected to slow 
down, as is the sunward drift of CAIs in the dead zone. The con- 
clusion is that CAIs can be transported to the dead zone before 
its formation. 



have departed from a sufficiently hot disk region, their transport 
range may be somewhat smaller than Rynax (e-g- ~10 AU in Fig. 
2). 



4. Preservation of CAIs in the dead zone 

In the self-similar model the return to the Sun, which es- 
sentially follows the gas, is completed for t - 4.9 fmax that 
is, largely before chondrite accretion. Although the estimates 
here are conservative, because they ignore random-walk turbu- 
lent motions (e.g. |Cuzzi et al.|2003[|Ciesla|2010| l, the problem is 
especially acute for centimeter-sized CAIs common in CV chon- 
drites (with = 10-100 kg/m^). The purpose of this section is 
to show that the appearance of a dead zone, which slows down 
the inward motion of the solids, solves the preservation issue. 

Indeed, thus far we have not addressed the question of the 
driver of accretion, and hence the maximum domain of valid- 
ity of the model. These drivers are thought to result from GI 
and/or the MRI. The red zone in Fig. 2 shows the maximal do- 
main (calculated ignoring heating by local turbulence) where GI- 
induced turbulence can occur, assuming a threshold value for the 
Toomre parameter Q = Qcs/nGl, of Qc - 2 (see e.g. 'Kratter| 
et al.|2008) l. Although GI may be important in the early phases, 
MRI-powered turbulence is believed to prevail after the disk has 
significantly expanded. The blue zone shows where neither tem- 
perature (below an activation threshold here fixed at 1500 K) nor 
cosmic rays (limited to an attenuation depth of lO"* kg/m-) yield 
the minimum ionization fraction required for the MRI (Gammie 



^ This result is obtained by numerically integrating the (Lagrangian) 
equation of motion of the solids at the last stagnation point, i.e. where 
and when the last positive velocities in the disk disappear. 



|1996[ ). Thus, wherever the blue zone is not overlaid by the red 
one, a values assumed by the self-similar model (10^^ - 10"^) 
cannot be self-consistently sustained, and a dead zone, with re- 
duced accretion, results within lO'* - 10^ yr at a few AU (stage 
'C in Fig. 1). This is before the CAIs have returned to the Sun. 

In the evolved stages of the disk (no longer captured by 
the self-similar model, because its simple a prescription breaks 
down), the inner edge of the dead zone will move inward to a 
fraction of an AU, while its outer edge will lie at tens of AU (e.g. 



Sano et al.|2000l|Ilgner & Nelson|2006[|Bai & Goodman|2009 1. 
The dead zone would then persist over most of the gas disk's 
lifetime: 'Fromang et al. (20021 and Terquem (2008) found e.g. 



dead zones in steady a-disks for a wide parameter range, in- 
cluding mass accretion rates at the lower end (~ 10^^ Mo/yr) of 
those of observed T Tauri objects ( Hartmann et al.|I998 1. 

Accretion is reduced in the dead zone, and the gas density is 
high enough to also reduce Vdrag (stage 'D' in Fig. 1). Using Eq. 
Q and (|6]l, the CAI drift timescale is 



fdrift = 



R 



V„ -I- Vdr; 



2MyrR 



1/2 
AU 



280 K 10-"* 
T max(Q', St) 
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The Stokes number St = npsajTL is related to the gas drag and a 
to advection. Note that this general order-of-magnitude estimate 
is independent of the previous history of the disk, a w ould not 
exceed 10 in the dead zone ( Fleming & Stone||2003 l and the 
surface density would likely be in the range 10 - 10^ kg/m^ 
( Terqueni||2008 1, so that millimeter-sized bodies satisfy St < 
0""*. This should ensure survival for a few Myr at a few AU. 
After the formation of the dead zone, the CAI production 
would essentially be confined to the inner active zone. Indeed, 
an important point is that the temperature of the thermal activa- 
tion of the MRI is on the order that of CAI condensation. Beyond 
the inner edge of the dead zone, the temperature will drop sig- 
nificantly below this value because of the drop of the turbulence 
level. Because diffusion from inside this boundary would be neg- 
ligible, CAIs formed after the dead zone appearance would not 
be sent to the chondrite-forming regions in significant amounts. 

5. Discussion 

Using viscous expansion of the disk to transport solids 
in general (Wehi-s tedt\fe Gail|[2002l [Gara ud '2007'. [Hughes &| 
Armitage 201 0) or specifically CAIs ( [Ciesla 2010) — with the 
recognition that their observed age distribution is narrowed be- 
cause CAI transport rapidly becomes inefficient — , is not new. 
However, the input parameters of the iTins of [Ciesla (2010 1 al- 
ways placed the CAI formation region inside /^switch and there- 
fore required diffusion to send CAIs to the spreading zone. In 
Sect. 2 we argued that an initial disk even more compact than 
envisioned by [Ciesla ( 2010 1 is conceivable and gave a prime role 
to advective transport. Its advantage is that the CAI-to-gas ratio 
would essentially retain its original order of magnitude, which 
seems required to account for the high mass fraction of CAIs 
in some chondrites. It would also explain the high proportion 
of crystalline silicates in comets ( Wooden et al.[[200F 



Hughes 



ow tur- 



& Armitage 2010). Finally, a small disk size would all( 
bulent mixing to achieve a significant initial degree of isotopic 
(and chemical) homogeneity in it, as is generally observed in 
meteorites (e.g. Cassen & Summers|1983 Yang & Ciesla 20T0) l. 
Regardless of the importance of advection, turbulent diffusion 
would help outward transport of CAIs and consequently, as far 
as the latter is concerned, the difference between this (very sim- 
plified) study and those of Bockelee-Morvan et al. ( 2002| l; Cuzzi 
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[eraL] ( |2003] l; [Ciesla| ( |20T0l ); |Hughes & Armitage| ( |20T0l i, is quan- 
titative, not qualitative, in nature. Because this model (or those 
listed just above) describes the origin and evolution of CAIs 
disk-wide, one prediction is that the key constraint of the re- 
stricted formation interval, hitherto drawn only from CV chon- 
drites, should pertain to CAIs of most if not all chondrite groups. 

We caution that infall should be taken into account in the 
future, although we ignored it in this study. This will require 
dedicated numerical simulations. Conceivably, outward expan- 
sion of the bulk of the disk beyond the centrifugal radius and 
appearance of the dead zone could be linked to a drop of the in- 
fall rate: matter supply would then be insufficient to sustain GI 
and would no longer overwhelm viscous-stress-driven flow past 
the centrifugal radius (e .g.|Zhu et al.||2010 In their protostel- 



lar collapse simulations, Tschamuter et al. ( 2009 [ ) also reported 



outward directed flows near the midplane in the hot inner disk, 
lasting a few millenia only (possibly because of their high as- 
sumed turbulence level), and mentioned the possible relevance 
for CAI formation. 

At the end of Sect. 3 we observed that the self-similar model 
fails to preserve CAIs when extrapolated over 1 Myr. This seems 
a general property of simple a-models. Ciesla (2010 1 found the 
discrepancy with observed CAI abundance to be between one 
and two orders of magnitude. He also noted that lower val- 
ues of a, while slowing down disk evolution, would also re- 



duce CAI production because of the lower temperatures. Cuzzi 
|et al.| ( |2003| l already noted this shortcoming and proposed that 
the inward drift of meter-sized migrators could maintain a 10- 
100-fold enhanced abundance of condensable matter above so- 
lar in the "CAI factory", hereby compensating for their subse- 
quent dilution in the nebula. However, this enhanced metallic- 
ity would yield partial pressures of oxygen in exce ss of those 
inferr ed from the Tv'^/Ti'^^ ratio in CAI pyroxene ("Grossmanl 
2010 i unless the migrators entered the CAI factory with C/O^ 
1 (see Eq. 24 in |Grossman et al.|2 008). This value lies far above 
that of the (chemically primitive) matrices of chondrites though 
and is consistent with the high end only of the range spanned by 



comet samples and interplanetary dust particles (e.g. Rietmeijer 



1998 J Sandford et al. 2008 ). Moreover, any enhancement would 



be limited by the amount of migrators in the outer disk ( Ciesla 
|2010 l), in particular for relatively compact disks. Ciesla ( 2010| l 



also argues that preferential settling of CAIs to the midplane 
could enhance their concentration in chondrites. With regard to 
the CAI drift timescale, it is unclear however when the low val- 
ues of the surface densities (~ 10^ kg/m^) he implies in his es- 
timate, which are two orders of magnitude below the MMSN at 
a few AU, would hold, especially in conjunction with low turbu- 
lence levels. 

Our point is that a solution to this problem naturally presents 
itself when the sources of the turbulence that drives mass ac- 
cretion are considered. Indeed, as pointed out in Sect. 4, simple 
a-models are not self-consistent throughout the disk extent and 
the evolution in the context of MRI- and/or GTdriven accretion. 
The formation of a dead zone is a general, unavoidable feature of 
accretion disks (Zhu et al. 2010| l: for instance, in a steady a-disk, 
the mass accretion rate (M = SttSv) below which Q right outside 
the inner MRI-active zone, evaluated with the passive irradiation 
temperature, exceeds Qc is given by 



M = 2xlO-M3/yr(— ) (-) ( 

^Ti „.(l AU) \''^V0.5 m'/kg 



1500 K 



5/8 



'Mo 



3/16 



280 K / 



(10) 



which is therefore a lower bound for the appearance of the dead 
zone in this case. Thus, the preservation mechanism outlined in 
Sect. 4 is actually fairly independent of the compact disk expan- 
sion scenario of Sect. 3. Ciesla (2010) mentioned among alter- 
native possibilities that a slowing of the disk evolution would 
explain CAI survival ; this study thus provides a way to realize 
this. That a dead zone blocks diffusion from the inner active zone 
would also account for the restricted formation interval of CAIs 
incorporated in chondrites, regardless of the dominating trans- 
port mechanism (advection, diffusion...) before the dead zone 
formation. The presence of a dead zone is also constrained by the 
Myr spread in age observed for chondmles (whose sizes are sim- 
ilar to that of CAIs) from individual chondrites (e.g. Villeneuvel 
et al.||20()9 1 by virtue of Eq. (j9|l. Finally, although we showed 
the necessity of the formation of the dead zone, how this for- 
mation occurs has yet to be ascertained. Because the viscosity 
drops outward around the inner boundary of the dead zone, a net 
outward flow may be expected there according to the first equal- 
ity in Eq. (|4| — perhaps even independently of the compactness 
of the disk — , before some steady state with uniformly inward 
accretion sets in. Detailed numerical simulations of dead zone 
formation are definitely needed. 
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